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Abstract. A treatise on the optical property of biological tissue is presented. Water is postulated 
to be a topological basis and serves to discriminate published skin data. Electromagnetic theory 
governing dielectric behavior is concisely detailed pertaining to certain optical constants and 
Kramers-Kronig relation. The Kramers-Kronig relation defining dispersion index is emulated 
through the discrete Hilbert transform. An accrued absorption set is populated with empirical 
absorption data for biological skin, pure liquid water and interpolated values. Kramers-Kronig 
analysis of biological skin yields a comprehensive description of the complex index of refraction 
from DC to x-ray frequencies. 
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1 INTRODUCTION 

A survey of published optical constants for biological skin underscores attention has been placed 
on the Radio Frequency (RE), near-infrared (NIR) and visible (VIS) frequency bands. Atten- 
tion has shifted recently to the terahertz band (THz), due to newly available radiation sources. 
Within each of these bands, reported values can differ dramatically from author to author and 
are commonly attributed to inherent error in experimental technique, numerical model used 
to calculate optical constants or random error. Experimentalists are limited to the narrow fre- 
quency band of their source and often cannot comprehend how their data fits into the broader 
picture. Theory can provide a comprehensive perspective of dielectric behavior, bring empirical 
data into focus and illuminate misconceptions. 

Biological skin is a multi-layered heterogeneous dielectric material. Three typical divisions 
of skin are epidermis, dermis and subdermis layers, where epidermis is frequently referenced as 
the stratum corneum. The partial molar configuration varies for each distinct layer of skin and 
also across racial groups, changes with age and is affected by disease, to name a few factors. 
Enumeration of so many variables discourages any unified treatment of biological skin. A global 
approach becomes accessible if water is postulated as a topological basis for all skin types. 
Consequently, it becomes possible to treat biological skin as an electromagnetic continuum. 

Electromagnetic theory describes the interaction of electric and magnetic fields with pon- 
derable media through constitutive relations. Each constitutive relation must be known for the 
bandwidth of radiation considered to accurately model and predict electrodynamic behavior of 
biological skin. Theory cannot produce results independent of experimental data. Electromag- 
netic theory permits calculating dispersion index from known absorption behavior through the 
well-known Kramers-Kronig relation j lj. Kramers-Kronig analysis of aggregate tissue absorp- 
tion data yields a comprehensive dispersion curve spanning the entire frequency spectrum. 
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2 THEORY AND METHODS 

Maxwell's equations, combined with all appropriate constitutive relations, fully describe elec- 
tromagnetic behavior in ponderable media 0][2]. There are three relevant constitutive relations 
for electromagnetic phenomena in dielectric media; the permittivity, admittance and perme- 
ability. Each relation is represented as a complex function of frequency and position, where 
dependence on position can be neglected in isotropic media. The complex permeability (/t) 
determines magnetization of a material in response to externally applied magnetic fields. The 
complex admittance (a) relates applied fields to induced current densities and is the sum of 
conductivity (er') and susceptivity (cr"). Conductivity measures the ease direct current flows 
through a material and susceptivity measures the ease dynamic current changes direction. The 
permittivity (e) measures the influence of an electric field to charge migration and dipole reori- 
entation. The real and imaginary parts of the complex permittivity describe stored energy in the 
system (e') and a dissipative term (e") measuring the loss of energy to the system. 

e = e' + ie" 

p, = n' + ifjf' (1) 
a = a' + ia" 

The complex index of refraction (N) combines all three constitutive relations into one 
complex function J2][3|. Because the complex index of refraction is a dimensionless func- 
tion, each constitutive relation is relative to vacuum, where the vacuum permittivity is eo ~ 
8.854 x 10~ 12 F/m and vacuum permeability is [Iq = An y. 10~ 7 H/m. For non-magnetic ma- 
terial, the relative permeability is considered equal to unity. The dispersion index («) measures 
frequency specific variation in phase velocity and accounts for reflection, refraction, chromatic 
dispersion and elastic scattering of incident electromagnetic radiation. Both Rayleigh and Mie 
scattering are examples of elastic processes causing redirection of electromagnetic energy. The 
absorption (k) term measures loss of energy to the system, including nonconservative losses 
through inelastic scattering processes. 



N(cj) = n + ik = \l fi r e r + i — r (2) 
V w 

Electromagnetic theory states absorption is null as frequency passes through the origin and 
infinite point Q]|2]- Equation (f3]i shows the one-sided limit for dispersion index and absorption 
in terms of elementary constitutive relations, where a = e' r — 4ira" /oj, b = e"+4:ira' r /u. Invok- 
ing Minkowski's inequality to approximate the inner square root term found in each definition 
reduces each one-sided limit, see Theorem 12. II The dispersion index is directly proportional 
to real permittivity, referenced as the static permittivity or static dielectric constant. Absorp- 
tion is shown primarily dependent upon the imaginary permittivity. This property of absorption 
enables an accurate calculation for low frequencies knowing only conductivity, where the fol- 
lowing relation exists: e" = a'/uj. 



lim < ^ = f + ^i U £->\ ) 



Theorem 2. 1 (Minkowski's inequality) The following inequality exists for two real 
numbers a and b ~ 

(a 2 + b 2 ) 1/2 < a + b, {a,be K} (4) 

This relation becomes an equality if either a or b equal zero. 

In the limit of u —> oo, the dispersion index approaches unity, while absorption rapidly 
approaches zero HIE). The physical reason for this behavior is that polarization has not enough 
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time to respond to rapidly changing fields E and the processes responsible for differentiation of 
the induction field D do not occur. The frequency where this behavior becomes valid is referred 
to as the plasma frequency (uj v ) and for light elements occurs somewhere in the far-ultraviolet 
spectra JTj. For frequencies higher than the plasma frequency, the dispersion index is solely 
dependent upon the concentration of conducting electrons, thus, the uniqueness of two similar 
materials become negligible. 

Corollary 2.2 (Non-uniqueness) The uniqueness of biological skin will fade in the 
limit of infinite frequency, hence, the absorption characteristics will approach that of pure liquid 
water. See also Postulate d. 6\ 

In general, the processes represented by term (a) tend to mutually cancel in the definition 
for absorption, but not in the case for dispersion index. Except for DC and infinite frequen- 
cies, the dispersion index is explicitly dependent upon all physical processes occurring in a 
material, including absorptive processes. Experimentalists employ various theoretical models 
to calculate a material's optical property based on measured reflectance and transmittance data. 
The total attenuation (/i t = fx s + fj, a ) of incident radiation through a material is identified as 
transmittance and is the combination of both the scattering coefficient (fi s ) and Lambert's ab- 
sorption coefficient (/i a ). From measured reflectance data and theoretical model employed, the 
magnitude of each coefficient is determined. This is a difficult proposition, especially for turbid 
media like biological tissue. For lower frequencies, effects such as electrode polarization and 
Maxwell- Wagner effect obscure measured results [0H6). For optical frequencies, experimental 
methods and theoretical models are often found inadequate irT HTOl . 

Absorption can be isolated experimentally, where the total attenuation can be simplified 
to Lambert's absorption coefficient (/i a ) alone. Since both absorption and dispersion index 
are spatially dependent, careful consideration of sample thickness and measured reduction in 
radiation intensity passing through the sample can be modeled with Beer-Lambert's law ifTTI . 
Lambert's absorption coefficient has a simple relationship to absorption (/i a = 47rfc/A), where 
A is the wavelength of radiation. The dispersion index shares no equally simple relationship in 
terms of fundamental constitutive relations. 

It is for the above reasons that some reliable method is needed to generate the index of re- 
fraction for biological skin. Fortunately, electromagnetic theory provides an analytic transform 
defining the dispersion index Q}0. The Kramers-Kronig relation defining the dispersion in- 
dex («) in terms of absorption (k) is shown in equation (0. This relation holds for any bound 
complex analytic function in the upper half plane. Analyticity in the upper half of the com- 
plex plane implies causality |fl][2][T2]- Formula is of particular importance, for it makes the 
calculation of the dispersion index, even approximately, from the absorption function. On the 
other hand, the uniqueness of the absorption function in terms of dispersion index is not guar- 
anteed, for the condition n(u) — 1 > is not satisfied in general. Therefore, implementation of 
Kramers-Kronig transform must be applied to absorption. 

n(u)-l = - r -^-du' = -— *k(u) (5) 

The limits of integration demand knowing the behavior of the absorption function for the 
entire real line. Causal systems are represented mathematically by Hermitian functions, which 
are comprised of an even real part and an odd imaginary part. Since the absorption function is 
bounded and odd, then by theory, the absorptive behavior of any material is null as the frequency 
passes through the origin and the infinite point. The absorption behavior for any material in the 
negative half-plane can be constructed by taking the negative of the positive half-plane and 
reversing the order along the abscissa. 

The convolution form of the Kramers-Kronig relation is shown in equation (0, where con- 
volution is signified by the asterisk symbol (*) lfT2l . In convolution form, the formula is more 
commonly referred to as the Hilbert transform. The Hilbert transform is most conveniently 
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evaluated with the aid of Fourier transforms. Applying the Fourier transform to a convolution 
results in the multiplication of the Fourier transform of each function separately. Converting the 
convolution integral in equation © to the transform domain requires taking the forward Fourier 
transform (J£") of the discrete absorption values for biological skin and then multiplying by the 
sign function (sgn), which is then multiplied by the imaginary number (i = \J — 1) lfl2l . The 
sign function is defined in equation (O and is the Fourier transform of the Hilbert transform 
kernel (— I/ttlo), where the transform takes cj n- s. To recover the dispersion index, unity is 
added after the inverse Fourier transform (J^" _1 ) is applied to the product. 

n(u) = 1 + sgn(s) J?{fc(u;)}} (6) 

The sign function is defined to be ~ 

C -1, s < 

sgn(s) = ^ 0, s = (7) 
[ 1, s>0 

The Hilbert transform (.ff) of a Gaussian function is shown in equation ([8]) and is a revealing 
identity. Since the absorption peak can generally be represented by a Gaussian function, the 
identity equals the real (5ft) part of the product between a Gaussian and error function (erf) 
with imaginary argument. The typical behavior for anomalous dispersion is controlled by the 
relaxation time constant (t), which also controls the width of the absorption peak. Both the 
prominence and steepness of the dispersion curve increases with increasing relaxation time. 

jr{exp(-rx 2 )} = 5R{i cxp(-rx 2 ) erf (iy/rx) } (8) 

Confusion may exist concerning the role relaxation time has with regard to frequency re- 
sponse and the shape of the anomalous dispersion curve. Consider the Similarity Theorem from 
Fourier theory, a dilation of the relaxation time corresponds to a contraction in the transform 
domain, see Theorem l2.3l [,12il. The relaxation time is dilated for bound water; consequently, a 
contraction of the thermal absorption peak is expected in the frequency domain for bound water, 
which leads to Corollary 12. 41 Conversely, the Hilbert transform has a one to one correspondence 
for the relaxation time, suggesting a more prominent anomalous dispersion centered about the 
thermal resonance for bound water, see Corollary 12. 5 1 

THEOREM 2.3 (SIMILARITY) The dilation of a constant r corresponds with a contraction 
in the Fourier transform domain, where t uj. 

&{f{Tt)}=F{w/T)/T (9) 

COROLLARY 2.4 (BOUND WATER) The relaxation time is dilated for bound water; conse- 
quently, a contraction of the thermal resonance peak for bound water and biological tissue is 
expected. 

COROLLARY 2.5 (ANOMALOUS DISPERSION) The relaxation time has one to one corre- 
spondence for the Hilbert transform; consequently, the anomalous dispersion is expected to be 
more prominent around the thermal resonance frequency for biological tissue. 

To ensure perfect reconstruction, Shannon's sampling theorem demands a band-limited 
function is sampled at least twice the rate of the highest frequency component in the signal IfOl . 
Otherwise, the signal will be under-sampled, which means fine structure represented by high 
frequency components will be poorly reconstructed. It is equivalent to applying a low pass fil- 
ter to a signal. Due to limitations in computer resources, under-sampling is often unavoidable; 
nevertheless, under-sampled data can produce acceptable results ifTZl . 

Since the discrete Fourier transform (DFT) will be utilized to evaluate the Hilbert transform, 
proper consideration must be given to the interval sampled and the placement of windows lfl2l . 
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The DFT is cyclical in nature and requires sampling the entire interval bounding a band-limited 
signal or a complete cycle or multiple of cycles for a periodic function; otherwise, boundary 
errors will be introduced into the calculation. The DFT also requires elements sampled at 
regular subintervals. Since experimental data is irregular, interpolation techniques must be 
employed to properly implement the Kramers-Kronig transform. 

Large intervals of missing empirical data present a more difficult problem. Short of produc- 
ing empirical data, there are three common methods used to attack large gaps in experimental 
data: 1) generating data points from a theoretical model, 2) supplanting absorption data from a 
similar material or 3) interpolation. No theoretical models exist for biological tissue enabling 
accurate absorption predictions for arbitrary frequency. Supplanting data from another mate- 
rial is prohibited for pure substances. In the case of heterogeneous materials, supplanting data 
is permissible for substances sharing similar chemical composition; therefore, would naturally 
share similar absorption characteristics. Lastly, interpolation is accurate and stable if the inter- 
val between known data points is small, but this method is problematic for larger intervals. 

Despite variation in concentration, water is present in all types of skin and represents a 
considerable percentage of biological tissue's composition. We postulate that biological tissue 
should share common absorption characteristics with pure liquid water, especially in the high 
frequency limit, see Corollary 12.21 Also, we recognize that tissue should diverge from the ideal 
absorption behavior of pure liquid water for very important theoretical reasons |[T1l51 fL2lfT4ll . An 
additional absorption peak should appear in the absorption spectra due to electrophoresis of ions 
in solution. Also, the absorption peak centered around the thermal resonance should contract, 
see Theorem l2.3l and Corollary 12. 41 The presence of melanin, hemoglobin and proteins should 
increase absorption in the visible and ultraviolet spectra fBl . 

Postulate 2.6 (Topological basis) Pure liquid water is an element of the basis set B 
defining biological skin. 

A basis B for topological space X generates a unique topology and generates alternate 
topologies when combined with other open sets lfl5l . The addition or variation of constituents 
in tissue is equivalent to a collection of sets comprised of the basis B set intersecting a series of 
open sets s 6 S, where S C X . New topologies are generated for each unique collection and 
maps all possible tissues, for every animal. 

The accuracy of interpolation decreases proportionately to the interval width IfTBI . Increas- 
ing the interpolation order (p) does not necessarily remedy the situation, as order p increases 
so does the implied order of derivatives increase. Inherent inaccuracies in data can lead to 
wildly oscillating results when attempting high order interpolations. Richardson's extrapola- 
tion algorithm takes a set of approximations and through a clever algebraic relation reduces 
the overall order of error to produce one accurate approximation |fl61 . For a set of known data 
points {xo, ■ ■ ■ ,Xi, . . . , x n }, Neville's iterated interpolation algorithm produces an approxima- 
tion x with implicit interpolation order p n IfTBI . The interpolation order is directly related to the 
number (n + 1) of known data points used to form the approximation. For N approximations 
generated with sequentially increasing interpolation orders {pi, . . . ,Pn}, we propose applying 
Richardson's extrapolation algorithm to produce a more accurate approximation. The exact 
source of error need not be specified to implement Richardson's algorithm, it is only required 
that the error be stepwise dependent, which we assume related to the interpolation order p. This 
implies an order N for Richardson's method. All algorithms were written and executed in Wol- 
fram's Mathematica software installed on an Intel® Core™2 CPU T700 @ 2 GHz computer 
running Microsoft® Windows XP™ with 2 GB of RAM. 

3 EXPERIMENT 

The Kramers-Kronig transform requires an absorption set A spanning the entire frequency 
spectrum be defined for biological skin. Surveying available experimental data for biologi- 
cal tissue reveals publication primarily in three frequency bands. Each experimental data point 
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was digitized from published graphs, see colored points in Fig. (Q]). Because of Corollary 
12.61 Segelstein's theoretical absorption curve for pure liquid water is included for immediate 
comparison |3). 




2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 

logio(Hz) 



Fig. 1 . Published absorption data for biological skin (colored points) and Segelstein's absorption 
curve for pure liquid water (solid black line). The width for each shaded region; Region I: DC 
to 22 GHz, Region II: 22 GHz to 3 THz, Region III: 3 THz to 2.5 fxm (IR-B), Region IV: 2.5 
fj.m (IR-B) to 189 nm (far-UV), Region V: 189 nm (far-UV) to x-ray. 

Experimental conductivity data, ranging from 20 Hz to 20 GHz, supplied by one definitive 
source populates Region I, see orange colored points in Fig. (Q]i [ 11). It was found necessary to 
multiply calculated absorption (k) values by the vacuum permittivity (eg) and angular frequency 
(cj) after converting reported conductivity values (cr' exp ) using equation (|2}. The correction fac- 
tor became evident after comparing calculated absorption values to that of pure liquid water. 
The additional correction factor is attributed to polarization of dielectric layers within skin sam- 
ples. The corrected relation for experimental conductivity removes an additional dependency 
on capacitance from the Maxwell-Wagner effect (eaujo-' exp = e" = <r' r /uj) HEED- Gabriel's 
wet skin data confirms Corollary 12.41 showing a contraction of the thermal resonance peak. An 
additional broad peak is seen centered about 6 log 10 (Hz), evidently from the influence of elec- 
trophoresis. Gabriel provides 171 data points with subinterval (1.2999, 10.3010) log 10 (Hz) to 
absorption set A. 

Because of the Similarity Theorem 12.31 and reasons of symmetry, a similar reduction in 
absorption is expected in the THz band. Of the three authors reporting absorption for biological 
skin, 38 data points from He et al. reported for rat skin are added to absorption set A, spanning 
subinterval (10.5235, 11.8942) log 10 (Hz) |fT9"l . Absorption values reported by Berry et al. for 
in vivo human skin spanned too small an interval and did not characterize the most important 
area centered around the thermal maximum 11201 . Pickwell et al. employed a double Debye 
permittivity model to calculate reported absorption values for in vivo human skin lETl . The use 
of a double Debye permittivity model to extract optical constants returned values very similar to 
pure liquid water. Debye's permittivity model neglects intermolecular interactions and cannot 
represent contraction of the absorption peak ll22l . 
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Fig. 2. Microwave-THz region: Compares reported absorption behavior of human skin (colored 
points) to pure liquid water (solid black line). Both Berry et al. and He et al. show a reduction in 
absorption for skin, while Pickwell et al. reported absorption values for human skin very close 
to pure liquid water. 

The presence of such scatterers as melanin, hemoglobin, and proteins are expected to raise 
absorption in the visible and ultraviolet spectra 1141 , Yet, the absorption behavior for biological 
tissue is expected to resemble water in the infrared spectra. To a lesser or greater degree, this 
expectation is met by published experimental values collected, see Fig. ([3). Unfortunately, 
large discrepancies amongst reported values demand further discussion to discriminate which 
are valid. 

Pickering et al. found the single integrating sphere (SIS) method can suffer error 60% higher 
than the double apparatus method 1101 . The double integrating sphere method was shown to 
return non-unique results and fails for highly scattering media. Absorption values reported by 
Hardy et al. show the presence of two vibrational peaks, but the entire curve appears flat and 
non-responsive |23l . Hardy et al. used a derivative form of the single integrating sphere method, 
goniometry, where a reflective sphere is replaced by a detector on a rotating arm. 

In addition to error associated with experimental technique are numerical models used to 
calculate optical properties from reflectance and transmittance data. Yang et al. determined the 
widely used Kubelka-Munk (K-M) model fails for low absorptive media J8). Authors using both 
the SIS technique and K-M method report the highest absorption values, specifically, Gemert- 
Anderson and Gemert-Wan, where data was digitized from graphs published by Gemert et al. 
l24l . Salomatina et al. reports lower absorption for ex vivo human skin, where optical constants 
were calculated using a Monte Carlo model in addition to using the SIS method [25 ] . 

Troy et al. employed the more accurate double integrating sphere apparatus and reports 
absorption values for ex vivo human skin in the infrared spectra similar to pure liquid water 
l26l . Arimoto et al, Bruulsema et al., Taroni et al. and Doornbos et al. used either time or 
spatially resolved reflectance spectroscopy methods on in vivo human skin samples 127 30l. Of 
particular interest are values reported by Simpson et al, where a unique application of the SIS 
method appears to remedy the inherent inaccuracies 1311 . Simpson et al. developed a unique 
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calibrating scheme to correct for systematic errors and reports values similar to Troy et al, 
Bruulsema et al, Taroni et al. and Doornbos et al. 

In light of Postulate d. 61 reported absorption for biological skin resembling pure liquid water 
cannot be accidental. All reported values agree skin's absorption diverges midway through the 
near-infrared spectra from pure liquid water, circa 14.5 log 10 (Hz). After considering insights 
from both Pickering et al. and Yang et al, reported absorption that is excessively high is judged 
invalid. The magnitude of skin's absorption in the visible spectra is on the order of -6 log 10 (k) 
and will contribute very little to the Kramers -Kronig integral; therefore, wrangling over reported 
values is not warranted. 94 points of data reported by Doornbos et al. are added to set A. This 
experimental set was chosen for two reasons, 1) the reported interval contains the point where 
skin's absorption diverges from pure liquid water and 2) the interval extends furthest into the 
visible spectra with range (14.4665, 14.7795) log 10 (Hz). 



wavelength - microns 
1.893 1.194 0.7536 0.4755 0.3 0.1893 




14.2 14.4 14.6 14.8 15 15.2 

logio(Hz) 

Fig. 3. IR-VIS-UV region: Compares absorption for human skin from eleven different authors 
(colored points) to Segelstein's water curve (solid black line). The use of SIS method and/or 
the Kubelka-Munk model leads to exceptionally high absorption. 



Because of similarities with pure liquid water and Postulate 12.61 it is redundant to average 
reported absorption for biological skin in most of the infrared spectra. The vast expanse of the 
mid-infrared spectra is completely devoid of any experimental data, see shaded Region III in 
Fig. (fl3. An intense literature search came up empty handed. In light of reduced absorption 
in the THz band, it is not exactly clear where to place the left endpoint for supplanted water. 
In other words, approaching the THz band (He et al) from the mid-infrared, skin's absorption 
must diverge at some point if it is to connect smoothly. 1741 data points from pure liquid 
water are added to set A with range (12.3198, 14.46063) log 10 (Hz). The left endpoint of the 
supplanted water is arbitrary and interpolation is expected to supply needed data within the 
inserted void. 

The high frequency limit behavior, Corollary 12.21 supports supplanting water absorption 
data for skin in the high frequency limit. To set A are added 20374 data points from water 
with range (15.2, 21) log 10 (Hz). It is not known how to connect the right edge reported by 
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Doorbos et al. with the left edge of the high frequency limit behavior. Both the single integrating 
sphere method and the Kubelka-Munk model were used to derive reported absorption in the 
ultraviolet spectra. Both methods have been determined inaccurate for low absorptive turbid 
media; consequently, another gap was left intentionally vacant in the experimental absorption 
set A with range [14.7795, 15.2] log 10 (Hz). 

Theory states material absorption is null at the origin. The amplitude is vanishingly small 
for frequencies left of Gabriel's data (u>/2ir < 20 Hz) and will not contribute greatly to the 
Kramers-Kronig integral. We join the origin to Gabriel's data using a simple linear regression 
model (y = ju>) that intercepts the origin, with calculated slope 7 = 1.93041 x 10 -8 . 1,000 data 
points are appended to set A with frequency range (-32, 1.2999] log 10 (Hz). The experimental 
absorption set A now possesses an overall range of (-32, 21) log 10 (Hz). Interpolation is required 
to both generate a regular set of data and remove two large voids, specifically, the first void is 
near the THz band, range [11.8942, 12.3198] log 10 (Hz), and the second within the ultraviolet 
spectra, range [14.7795, 15.2] log 10 (Hz). 

From set A a new set A' is created using a combination of interpolation followed by ex- 
trapolation. Approximations were made using Neville's interpolation algorithm for arbitrary 
frequency ojq. The interpolation order p swept through all even integers 2n to produce a total 
of 30 approximations, where n = {1, 2, . . . , 30}. Richardson's extrapolation algorithm, with 
order 2n = 60, was then successively applied to all 30 approximations to yield a single accurate 
approximation x, for frequency ojq. This process was repeated for all 10,000 chosen frequencies 
regularly spanning the experimental spectrum (-32, 21) log 10 (Hz). Computational time was just 
under 9 hours to evaluate all 10,000 elements to form set A'. The algorithm was applied in base 
ten logarithmic space to minimize numerical error. 
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Fig. 4. Interpolated/extrapolated absorption set A' for biological skin (solid black line) and 
Segelstein's water absorption curve (dashed line). Inset focuses on the NIR-VIS-UV spectra. 
For details on lettered arrows see text below. 

Fig. © shows the resulting absorption curve A' for biological skin after interpolation and 
extrapolation. An experimental artifact is retained from Gabriel's data, see arrow A in Fig. ©. 
An extra peak is introduced by the interpolation scheme over the first void (arrow D). Inter- 
polation filled the void intentionally left in the ultraviolet spectra, where the right endpoint of 
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Doornbos et al. (arrow B) connects to the left endpoint of supplanted water values (arrow C). 
To test any possible ill effects in forming set A', measured relaxation times before and after in- 
terpolation/extrapolation are compared. Measuring the Full Width Half Maximum (FWHM) of 
the thermal peak from both absorption sets, A and A', returned the same approximate relaxation 
time, specifically 16.824 ps for skin. For comparison purposes, we measured the relaxation time 
from Segelstein's pure liquid water curve as 5.38 ps. 

The Hilbert transform was used to emulate the Kramers-Kronig relation, see equation (0. 
To satisfy Shannon's sampling theorem for set A', a sampling rate of 5 x 10~ 22 is required 
which would create a set with 2 x 10 21 point density. This far exceeds available computing 
power and forced implementing two under-sampled windows. The absorption magnitude is 
vanishingly small in the visible spectrum and provides a perfect place for the right endpoint of 
the first window, range (0, 14.8) log 10 (Hz). It is satisfactory to place the edge of the second 
window far into the x-ray band where absorption is vanishingly small, range (0, 17) log 10 (Hz). 
The windows are judiciously spaced to minimize boundary errors from the discrete Fourier 
transform (DFT). 

Based on available computing power, a point density of 2 21 was chosen for each window 
submitted. Proper implementation of the DFT requires a regularly sampled set of data points 
and each window shifted the placement and spacing of the resulting subintervals. This requires 
another round of interpolations to be performed on set A'. The spacing between each data 
point was considered sufficiently close to allow first order interpolation. Interpolations were 
performed in base ten logarithmic space utilizing Mathematica's built-in Interpolation!] algo- 
rithm [32 1 . The first approximation is always taken one step-size from the origin. Interpolation 
produces a list k with length 2 21 for each window, where the subscript i = {1,2} identifies 
which window is referenced. In order to represent absorption for negative frequencies, a new 
list l\ is created by taking the negative of list U in reversed order. Combining lists U and inverse 
lists l[ creates a composite list d representing skin absorption for both positive and negative 
frequencies. Below is a concise definition of the composite list Ci using Mathematica's notation 
and functions: 

d = Flatten[{-Rest[Reverse[ k ]], 0, k}] (10) 

Theory states absorption is null at the origin, hence, the number zero is explicitly stated. 
The function Flattenf] removes all nested lists within the outer curly brackets. The function 
Rest[] removes one element from the list enclosed and this is to maintain an overall length for 
the composite list Ci equal to 2 22 . Lastly, the function Reverse[] obviously reverses the order 
of elements in the enclosed list. 

Using the definition for the sign function (sgn), another list £ sgn of length 2 22 is created 
spanning both positive and negative frequencies, see equation @. Using Mathematica's built-in 
Fourierf] function, the forward discrete Fourier transform ( j£") is taken of the composite list d 
and then multiplied by list £ s i gn . The discrete inverse Fourier transform (i^ -1 ) is taken of the 
product list, ,^ r ~ 1 {I ■ £ sgn x J?{d}}, using Mathematica's built-in InverseFourierf] function, 
where I is the imaginary number V - 1- Mathematica's discrete Fourier transform functions 
will implement Cooley and Tukey's Fast Fourier Transform (FFT) algorithm for lists having 
length equal to an integer power of 2 11331 . The dispersion index is retrieved by adding unity to 
the list after the inverse Fourier transform is taken. Runtime was approximately 20 seconds to 
complete all numerical operations and process the 4,194,304 elements for each window. Half 
of the dispersion list generated is discarded, for we do not need the negative frequencies. 

Segelstein emphasized in his thesis that calculated dispersion curves are shape preserving, 
but not their amplitudes 0. This is a consequence of under-sampling, which poorly represents 
the area under the absorption curve. The amplitude is commonly corrected by scaling calculated 
dispersion to an excepted static permittivity value for the material under study. If a series of 
windows are used, then the disperison index from the first window is scaled and each subsequent 
output is scaled to match the right endpoint of each preceding corrected output. 
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In the present study, the static permittivity is not known for biological skin. In the limit 
of static fields, the dispersion index is directly proportional to the square root of the real per- 
mittivity (n 2 = e' r ), see equation (f3]). It is common to see values with orders of magnitude 4, 
5 and as high as 6 reported for skin's real permittivity l5ll7l [T7l[34l ■ An estimate of the com- 
pound dielectric constant for heterogenous mixtures is attainable through the mixing rules and 
Maxwell-Garnett's formula ll35ll . Pethig et al. reports static permittivity values of 90 and 110 
for amino acids and proteins, respectively 0. Since water represents roughly 75 percent of 
tissue, fractional percentages for amino acids and proteins of 10 and 15 percent return a value 
of 83 for the effective permittivity. Gabriel reports a real permittivity (e' r ) of 56432 for wet 
skin at 20 Hz, this stands in stark contrast to estimates obtained through the mixing rules and 
Maxwell-Garnett's formula. 

Because of the confusion concerning skin's static permittivity, another means of scaling 
dispersion index calculations is needed. It was noticed that predicted index values obtained 
from the second windowed transform compared favorably with index values reported by Ding 
et al. 11361 . Predicted values were within reported standard deviations for skin epidermis and 
dermis; moreover, the mean value approached very close to calculated index values. In light of 
this discovery, it was decided to reverse the order of scaling by scaling the output from the first 
window to match an arbitrary point in the second windowed transform. 

Wavelength 600 nm was arbitrarily chosen as the fixed point for scaling, the index predicted 
from the second windowed transform is approximately 1.41. This required multiplying the first 
windowed output by approximately 1.41 to raise the predicted value for 600 nm from unity. 
It is known that windowed transforms slope downward towards higher frequencies as a result 
of under-sampling the Hilbert kernel, it was chosen to multiply rather than add the correction 
factor for this reason lfT2ll . The resulting dispersion curve is shown together with pure liquid 
water's theoretical dispersion curve for reference, see Fig. Q. 




8 9 10 11 12 13 14 15 16 17 

logio(Hz) 

Fig. 5. Thoeretical dispersion index curve for both biological skin (solid line) and Segelstein's 
pure liquid water (dashed line). Both curves were generated through the Kramers-Kronig rela- 
tion defining index as a function of absorption. Inset focuses upon IR-VIS-UV spectras. 

By inspection, the dispersion index for biological skin is not strictly greater than pure liquid 
water. The prominence of the dispersion curve around the thermal resonance verifies Corollary 
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12.51 also, dilation of the relaxation time for bound water momentarily forces dispersion to dip 
slightly below that of pure liquid water, circa 10.8 log 10 (Hz). The predicted static permittivity 
is n 2 w 80.5 and compares well with estimates from Maxwell-Garnett's formula. The exper- 
imental absorption values collected are derived primarily from skin tissue samples, therefore 
the predicted dispersion is indicative of the combined index for the epidermis and dermis. Pre- 
dicted index values compare very well with the mean index obtained from averaging dermis and 
epidermis index reported by Ding et al. Il36l . 

4 CONCLUSION 

Kramers-Kronig analysis of the optical properties for biological skin produced a comprehen- 
sive set for the complex index of refraction spanning the practicable frequency spectrum, see 
Tables ((Tj-3). Segelstein's theoretical curve for pure liquid water proved an invaluable resource 
in discriminating experimental absorption for biological skin. Through comparison, an addi- 
tional correction factor was found necessary for experimental RF data; otherwise, the predicted 
slope was opposite to expected absorption behavior, increasing towards the origin instead of 
decreasing. Biological skin exhibits strong similarity to pure liquid water in the near-infrared 
spectra. Through comparison, derived absorption from the single integrating sphere method 
and Kubelka-Munk model were identified as exceptionally high and were not included in this 
analysis. 

It is surprising to discover results derived from widely accepted methods as the single inte- 
grating sphere method and Kubelka-Munk model were inaccurate. This should focus attention 
on the importance of global analysis, as well as, the need for further skin research in the mid- 
infrared and ultraviolet spectra. Of particular interest are the results reported by Simpson et 
al, the unique experimental method he describes for the single integrating sphere method ap- 
parently makes viable what is otherwise an inaccurate method when applied to low absorptive 
turbid media. 

A novel application of Richardson's extrapolation algorithm theoretically allowed greater 
accuracy for interpolated approximations. The artifact peak covering the intentional void near 
the THz band indicates a possible discontinuity in the tangent space between experimental skin 
data and supplanted water data. Due to a lack of experimental data, the exact behavior of biolog- 
ical skin is not known in the mid-infrared band. It is quite possible that we were too aggressive 
extending the left endpoint for supplanted water. In the case of the ultraviolet band, the interpo- 
lation/extrapolation algorithm produced a smooth regular curve connecting experimental data 
in the visible spectra with supplanted water used to represent the high frequency limit behavior. 
The implementation of the Kramers-Kronig transform requires only a simply connected curve. 
The absorption magnitude in the visible spectra is quite small and contributes little to the sum 
evaluated from the Kramers-Kronig integral. Even though the ultraviolet spectra is purely a 
product of interpolation, small changes in absorption values would not drastically change pre- 
dictions found in this analysis. The area in the ultraviolet spectra is not as important as the 
area around the thermal maximum. The thermal resonance peak represents the area of largest 
contribution to the integral and it is fortunate to find experimental values for this region. 

Comparison of biological tissue with pure liquid water is based on the assumption water 
is a topological basis set B. Postulate 12.61 imparts a global perspective and provides a more 
secure footing for analysis. With this concept, all possible biological tissue types are new 
topologies X created from the intersection of the basis set with other open sets (s 6 S). Open 
sets containing elements such as proteins, amino acids, melanin and others give rise to new 
topologies. Ultimately, there are a family of curves based upon partial molar concentration, 
specific elements contained in open sets and the degree of intersection with the basis set. The 
manifold M. of all possible tissue types is covered in a continuous fashion by the basis set and 
all included open sets. 
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To fully characterize all possible tissue types or even biological skin alone would require 
a comprehensive analysis including the variation of all pertinent parameters. This would be a 
daunting task and it may be more efficacious to treat the manifold Ai as a probability space. 
Segelstein's theoretical curve brought all experimental data into proper focus. In like manner, it 
is hoped the published theoretical curve for biological skin will aid future investigations. 
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Table 1 . The complex index of refraction for biological skin. 
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Table 2. The complex index of refraction for biological skin. 
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Table 3. The complex index of refraction for biological skin. 
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Table 4. The complex index of refraction for biological skin. 
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